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Abstract 

Background: The specific deposition of histone variants into chromatin is an important epigenetic mechanism that 
contributes to gene regulation through chromatin architectural changes. The histone variant H2A.Z is essential in 
higher eukaryotes, and its incorporation within chromatin is a relevant process for gene expression and genome 
stability. However, the dual positive and negative roles of H2A.Z in gene regulation still remain unclear. We 
previously reported that acclimatization in common carp fish {Cyprinus carpio) involves cyclical seasonal gene 
reprogramming as an adaptation response to its natural environment, when rRNA synthesis and processing are 
profoundly affected. Epigenetic mechanisms primarily contribute to the transcriptional modulation of ribosomal 
genes concomitant with the acclimatization process, thus significantly regulating this process. The aim of this study 
was to describe the presence of several H2A.Z subtypes in carp, and assess the role of H2A.Z on the ribosomal 
cistron in summer- and winter-acclimatized carp. 

Results: This paper reports for the first time about the transcriptional expression of four different H2A.Z subtypes 
belonging to the same organism. Remarkably, a novel H2A.Z.7 was found, which corresponds to a tissue-specific 
histone subtype that contains seven amino acid residues longer than the canonical H2A.Z. Moreover, H2A.Z 
enrichment through the ribosomal cistron was significantly higher during summer, when rRNA transcription and 
processing are highly active, than it was in winter. Similar patterns of H2A.Z enrichment are found in two seasonally 
active promoters for genes transcribed by RNA polymerase II, the L41 and A 9 -desaturase genes. Interestingly, 
ubiquitylated-H2A.Z (H2A.Zub) was strongly enriched on regulatory regions of the ribosomal cistron in summer- 
acclimatized carp. Additionally, H2A.Z was present in both heterochromatin and euchromatin states on ribosomal 
cistron and RNA polymerase II promoters. 

Conclusions: Our study revealed seasonally-dependent H2A.Z enrichment for active ribosomal cistron and RNA 
polymerase II promoters during the carp environmental adaptation. Moreover, seasonal H2A.Zub enrichment appears as 
a specific mechanism contributing to the regulation of chromatin architecture under natural conditions. The existence 
of several H2A.Z subtypes in carp suggests that the epigenetic regulation in this species constitutes a complex and 
finely tuned mechanism developed to cope with seasonal environmental changes that occur in its habitat. 
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Background 

Living organisms are constantly exposed to environmen- 
tal fluctuations. In order to survive, they must adapt by 
inducing compensatory responses [1]. In particular, tele- 
ost fish are an interesting model for studying environ- 
mental plasticity due to the fact that they are constantly 
exposed to several types of stressors in the environment. 
Thus, fish are living under large fluctuations of physical 
parameters, such as photoperiod, water temperature, nu- 
trition availability, oxygen concentration, and salinity, 
among others [2]. In this context, the common carp pro- 
vides a powerful, natural, non-manipulated model for 
studying several molecular and cellular responses as part 
of the acclimatization process [2,3]. 

It is well known that eukaryotic gene regulation and ex- 
pression are connected to the structural organization of 
chromatin. Chromatin, through its capability to transiently 
pack and unpack the genome, plays a key role in control- 
ling access to genetic information [4]. The basic unit of the 
chromatin is the nucleosome, which is formed by two cop- 
ies of each histone protein (H2A, H2B, H3, and H4) as- 
sembled into an octamer that has 145 to 147 base pairs 
(bp) of DNA wrapped around it to form a nucleosome core 
[5]. In addition to packaging DNA, chromatin is a dynamic 
structure that adjusts to different cellular processes in re- 
sponse to biotic and abiotic environmental parameters [6]. 
This flexible and efficient intercommunication between 
genomic DNA and exogenous influences can be explained 
by several epigenetic mechanisms, such as DNA methyla- 
tion [7], ATP-dependent chromatin remodeling [8], non- 
coding RNAs [9], post-translational modifications (PTMs) 
of histones, and the replacement of canonical histones with 
non-allelic histone variants [10,11]. Among these, the non- 
random deposition of histone variants into chromatin that 
occurs independently of DNA replication, substantially 
contributes to gene regulation and architectural changes in 
chromosomes [12]. 

H2A is the canonical histone that displays a great num- 
ber of variants, such as macroH2A, H2A.Bbd, H2A.X, and 
H2A.Z [13]. Among these, H2A.Z is the most evolutionar- 
ily ancient one between species, since it is highly con- 
served during evolution and is widely distributed in all 
eukaryotic lineages [14], The identity of H2A.Z compared 
with the canonical histone H2A is around 60%, suggesting 
a differential and particular role of this histone variant 
[15,16]. Indeed, H2A.Z genes are essential for the develop- 
ment and viability of many organisms, such as Drosophila 
melanogaster [17], Tetrahymena thermophila [18], Xenopus 
laevis [19], and mice [20,21], Additionally, absence of H2A. 
Z in Saccharomyces cerevisiae leads to slow growth and a 
deficiency in transcriptional induction by RNA polymerase 
II [22], In vertebrates, two non-allelic genes exist (with a 
common origin in early chordate evolution) that code for 
two highly similar H2A.Z proteins, termed H2 A.Z.I and 



H2A.Z.2 [23]. Even more recently in humans, two alterna- 
tively spliced variants for H2A.Z.2, named Z.2.1 and Z.2.2, 
have been described [24]. 

Several PTMs have been reported for H2A.Z, such as 
acetylation [25], sumoylation [26], and ubiquitylation [27]. 
In particular, the ubiquitylated form of H2A.Z (H2A.Zub) 
is predominantly found in the inactive X chromosome of 
female cells, which suggests that this PTM is involved in 
the formation of silent chromatin. Therefore, the factors 
that control ubiquitylation of H2A.Z, such as RING lb E3 
ligase [28] and Bmil [29], as well as those that regulate 
deubiquitylation, such as USP10 [30], have arisen as im- 
portant regulators for the epigenetic control of H2A.Z. 

rRNA transcription is a demanding and crucial mechan- 
ism for all eukaryotic cells because growing cells require a 
continuous ribosome supply [31]. Particularly in carp, this 
process undergoes dramatic changes in the context of 
adaptation to summer and winter [3]. Previously, we dem- 
onstrated that epigenetic mechanisms contribute to the 
transcriptional regulation of ribosomal genes concomitant 
with the acclimatization process. In this context, we have 
reported the hypermethylated state of ribosomal gene pro- 
moter in winter-acclimatized carp, as well the histone H2A 
replacement by variant subtypes of macroH2A, help the 
seasonal transcriptional activity of ribosomal genes [32,33]. 
In addition to trans-acting factors, we also reported on the 
genomic organization of the carp ribosomal cistron [34], 
when cis-acting elements located in the intergenic spacer 
(IGS), named T 0 ' and T 0 , and core promoter (CP), can 
contribute to the modulation of rRNA synthesis. 

Genome-wide location analyses have reported a non- 
random distribution of histone H2A.Z. However, despite 
intense efforts of many groups, H2A.Z function remains 
enigmatic. Previously, it was shown that H2A.Z is glo- 
bally localized in promoter regions [35]. In human and 
Drosophila cells, the enrichment of H2A.Z at gene pro- 
moters is positively correlated with transcription [36,37], 
whereas in yeast cells H2A.Z content of the same re- 
gions is generally inversely correlated with transcription 
levels [38,39]. 

Taking the above into account, and the limited evi- 
dence concerning the epigenetic response as a whole 
with a model under natural conditions, we evaluated the 
role of H2A.Z on the ribosomal cistron in summer- and 
winter-acclimatized carp. We found significant changes 
of H2A.Z deposition on ribosomal cistron depending on 
the seasonal transcriptional activity of rRNA genes. Fur- 
thermore, because PTMs of H2A.Z are associated with 
different chromatin states [27,40,41], we also found a 
differential enrichment of H2A.Zub in the IGS and CP 
regions. In addition, we report for the first time a tissue- 
specific subtype of H2A.Z, called H2A.Z.7, which is 
characterized by encoding seven additional amino acid 
residues, more than canonical H2A.Z. 
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Results 

Seasonal temperature and photoperiod variation in the 
natural habitat of the common carp in Santiago, Chile 

In order to determine periods when abiotic conditions were 
extreme, we analyzed the temperature and photoperiod as 
parameters that continuously fluctuate during the year 
2011. We found that the months between November and 
February (summer period) corresponded to the months 
with longer daylight (13 to 14 hours), and higher mean 
temperatures (17 to 22°C). Conversely, the period between 
May and August (winter period) corresponded with shorter 
daylight (10 to 11 hours), and lower mean temperatures 
(5 to 7°C) (Additional file 1: Figure SI). 

Isolation of different H2A.Z subtypes, and evaluation of 
tissue-specific expression of H2A.Z and H2A.Z.7 subtypes 

We isolated four different sequences that encode for 
different subtypes of H2A.Z in carp (Additional file 2: 
Figure S2). The deduced amino acid sequences were 
obtained and then compared with canonical H2A.Z se- 
quences previously described in the GenBank (Figure 1). 
The first carp sequence corresponded to H2A.Z.2.1 (100% 
identity with human H2A.Z.2.1). The second and third se- 
quences found corresponded to two new subtypes, which 
we named H2A.Z.3.1 and H2A.Z.3.2. Both subtypes differed 
in only one residue in position 11 (G and S, respectively). 
The fourth sequence was longer than H2A.Z and contained 
seven additional amino acids (positions 66 to 72). We 
named this subtype H2A.Z.7 according to new histone 
variant nomenclature [42]. 

With the purpose of evaluating tissue-specific gene ex- 
pression of different H2A.Z subtypes in seasonally adapted 
carp, RT-PCR assays were performed on different tissues. 
As shown in Figure 2A,B, the presence of the H2A.Z.7 
transcript was detected only in the brain and pituitary tis- 
sues. Particularly, in the brain, no significant differences in 
the transcriptional expression of H2A.Z. 7 were observed 
between summer and winter seasons. In contrast, the pitu- 
itary H2A.Z.7 transcript was only detected during the 



winter. On the other hand, we performed RT-PCR assays 
with primers that did not discriminate between H2A.Z 
subtypes (H2A.Z2.1, H2A.Z.3.1, and H2A.Z3.2). This 
amplification was denoted as H2A.Z, and it was detected 
in all examined tissues (Figure 2B). Higher levels of the 
H2A.Z transcript were detected in testis of summer- 
acclimatized fish, while the same tendency was observed 
in the liver during winter- acclimatized fish. No significant 
variation in the content of H2A.Z was observed between 
seasons in the brain, pituitary, and gills. 

The phylogenetic relationship of carp H2A.Z subtypes 

In order to analyze the phylogenetic relationship among 
the new H2A.Z subtypes and the previously described 
H2A.Z, we constructed a phylogenetic tree using 29 cod- 
ing DNA sequences (CDS) from the GenBank, including 
H2A.Z from three different kingdoms, such as Animalia, 
Fungi, and Protist (Figure 3). According to the phylogen- 
etic reconstruction, subtypes of the histone variant H2A.Z 
isolated from carp, together with the group comprising 
H2A.Z variants belonging to the Actinopterygii class of 
fish, share a monophyletic origin. Furthermore, these se- 
quences cluster together with the H2A.Z from other tele- 
ost fish, such as Danio rerio, Salmo salar, Pagrus major, 
Oncorhynchus rnykiss, and Esox lucius. 

The analysis revealed that H2A.Z subtypes from carp, 
and especially H2A.Z.7, are subject to a marked diversifi- 
cation process within this vertebrate group of animals. 

Nuclear content of the histone H2A.Z varies seasonally in 
carp hepatocytes 

In order to validate and confirm that commercial anti- 
bodies developed against mouse H2A.Z cross-react with 
samples of carp hepatocytes, comparative western blot 
with rat and human cells nuclear extracts was performed. 
We found a specific signal of approximately 14 kDa corre- 
sponding to carp H2A.Z, indicating that antibodies detect 
their orthologs in carp (Figure 4A). 
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Figure 1 Multiple amino acid sequence alignments between H2A.Z subtypes isolated from carp and H2A.Z variants from human. 

Protein alignment was performed using ClustalW software (Conway Institute, University College Dublin (UCD), Dublin, Ireland). Amino acid 
differences are highlighted with blue background. The additional residues of variant H2A.Z.7 isolated from carp are depicted in yellow. Protein 
sequences from human H2A.Z subtypes were extracted from Bonisch et al. [24]. 
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Figure 2 Tissue-specific transcriptional expression of H2A.Z and H2A.Z.7. (A) RT-PCR analysis of H2A.Z.7, H2A.Z and 16S mitochondrial rRNA 
(16S) from brain, pituitary, liver, testis and gills tissues from summer- and winter-acclimatized carp. Three different fish for each season were 
examined (S 1; S 2/ S 3 and W 1f W 2 , W 3 ). The last line (-) corresponds to the negative control of PCR reaction. (B) The graph shows a relative 
quantification of tissue-specific transcriptional expression of H2A.Z subtypes. Orange bars denote H2A.Z.7 and black bars represent the content of 
the remaining H2A.Z subtypes, respectively. All bars represent the ratio between the respective mRNA and carp 16S mitochondrial rRNA as a 
control. S, summer; W, winter. 



Subsequently, we analyzed the protein content of 
H2A.Z in hepatocytes from seasonally adapted carp. We 
found that in winter, H2A.Z is 1.7 times higher com- 
pared to summer-adapted carp (Figure 4B,C), a result 
that is consistent with the liver tissue-specific RT-PCR 
assay shown in Figure 2. 

Since the extracellular matrix contributes to the autoflu- 
orescence emission in hepatocyte tissue, we performed an 
H2A.Z immunodetection assay in pituitary tissue and 
found similar results as those previously observed in liver. 
Unexpectedly, nuclear protein content of H2A.Z from pi- 
tuitary cells was much higher in winter- acclimatized fish 
(Figure 4D). 

Seasonally-dependent distribution of H2A.Z throughout 
the ribosomal cistron in the carp 

In order to evaluate the possible differences of H2A.Z 
enrichment along the carp ribosomal cistron during 
acclimatization, we used chromatin immunoprecipita- 
tion (ChIP) assays and subsequent quantification by 



quantitative PCR (qPCR) (Figure 5A). The results showed 
a sevenfold increase (± SD) of H2A.Z enrichment in IGS, 
T 0 ', T 0 , and CP regions, and a threefold increase (± SD) in 
the 5.8S region during summer (Figure 5B). A H3 ChIP 
control to determine nucleosome occupancy during sea- 
sons at investigated loci is shown in Figure 5D. 

Additionally, with the purpose of correlating the differ- 
ential seasonal enrichments of H2A.Z with general gene 
transcriptional activity, we evaluated the association of 
H2A.Z with the promoter regions of three RNA polymer- 
ase II-related genes. Two of them, L41 and A 9 -desaturase 
genes, exhibit seasonally-dependent expression in carp, 
whereas JAK1 gene has a constitutive expression between 
seasons (see Additional file 3: Figure S3) (Figure 5C). Pre- 
viously, it has been described that L41 and A 9 -desaturase 
genes display increased transcriptional activity in hepato- 
cytes from summer- and winter-acclimatized carp, re- 
spectively [43,44], In agreement with this observation, 
H2A.Z was systematically enriched (by approximately 
three times) on promoters when genes are seasonally 
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Figure 3 Phylogenetic relationships among the histone H2A.Z family from different eukaryotes. The phylogeny reconstruction was carried 
out using the minimum evolution (ME) method and the Close-Neighbor-Interchange (CNI) algorithm. The root of the tree is denoted in black and 
the numbers represent p-distance (see methods for details). The four new carp H2A.Z subtypes are indicated in red. The sequences of H2A.Z are 
listed in Additional file 5: Table S2. CNI, Close-Neighbor-Interchange; ME, minimum evolution. 



activated. In contrast, on the promoter of the JAK1 -consti- 
tutive gene, there is no significant difference of H2A.Z en- 
richment between seasons. 

Deposition of H2A.Zub in summer-acclimatized carp 
throughout the ribosomal cistron 

The monoubiquitylation of H2A.Z distinguishes its asso- 
ciation with euchromatin or facultative heterochromatin 
[27]. Consequently, we decided to study whether or not 
there is a relationship between the H2A.Zub deposition 
on ribosomal cistron with the transcriptional activity of 
carp rRNA. The results showed that during summer, 
there was an increase of fourfold, threefold and sixfold 



(± SD) of H2A.Zub enrichment in IGS, T 0 and CP regions, 
respectively (Figure 6A); whereas, H2A.Zub enrichment in 
5.8S was not significant. To determine whether these find- 
ings are exclusive in ribosomal cistron or not, we evalu- 
ated the association of H2A.Zub with the promoter 
regions of the ribosomal protein L41 and A 9 -desaturase 
genes. In both promoter regions, H2A.Zub enrichment 
was not significant (Figure 6B). 

Moreover, we performed RT-qPCR assays against Bmil, 
RinglB, and USP10 ubiquitylation/deubiquitylation factors 
to reveal any possible variations in the ubiquitylation ma- 
chinery concurrent to seasonal changes (Figure 6C). The 
mRNA levels of Bmil and RING lb were fourfold and 
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Figure 4 Nuclear content of H2A.Z from hepatocytes and pituitary cells. (A) The cross-reactivity of commercially available antibodies 
between H2A.Z from rat, carp, and human cells (30, 100, and 50 ug per line of nuclear extract, respectively) is shown. (B) Western blot assays of 
H2A.Z and H2B in hepatocytes isolated from four different fish for each season (S 1; S 2 , S 3 , S 4 and W 1; W 2 , W 3 , W 4 ). (C) Seasonally relative 
quantification of H2A.Z content. The H2A.Z content was normalized against H2B content (n = 4 for each season). Standard deviations (± SD) are 
shown (*P = 0.340, Student's f-test). (D) Seasonal immunodetection of H2A.Z is shown from pituitary sections: (A to D), DNA detection with DAPI; 
(B to E), H2A.Z immunodetection; and (C to F), merge. Bar = 25 urn. DAPI, 4',6-diamidino-2-phenylindole; S, summer; SD, standard deviation; 
W, winter. 



sixfold higher, respectively, in the summer compared to 
winter. Conversely, there were no changes in the mRNA 
content of USP10 between seasons. 

Correlation between the enrichment of H2A.Z and 
epigenetics marks 

We define colocalization as the presence in the same 
chromatin region of H2A.Z, and either euchromatin 
(acetylated histone H4 at lysine 12 (H4K12ac)) or het- 
erochromatin (trimethylated histone H3 at lysine 9 
(H3K9me3)) markers. Using this colocalization, we can 
correlate H2A.Z occupancy with chromatin transcrip- 
tional states of the rDNA cistron and two RNA polymer- 
ase II promoters. 

We found a similar proportion of colocalization be- 
tween H2A.Z and H4K12ac between seasons along the 
rDNA cistron (Figure 7A). In contrast, when we ana- 
lyzed the colocalization between H2A.Z and H3K9me3, 
this was significantly greater during winter on the T 0 ', 
T 0 , and CP regions (Figure 7B). 

H2A.Z enrichment was higher on the promoter of the 
A 9 -desaturase gene during winter, co-localizing with the 
euchromatin marker H4K12ac. The same colocalization 



between H2A.Z and H4K12ac was found in the summer- 
active L41 promoter (Figure 7C). In contrast, H2A.Z en- 
richment co-localized with the heterochromatin marker 
(H3K9me3) when the transcriptional activity of gene ex- 
pression is seasonally repressed (Figure 7D), as in the case 
of L41 during winter and A 9 -desaturase during summer. 

Discussion 

In carp, seasonal acclimatization involves a reprogram- 
ming of molecular and cellular functions, which begin 
by sensing environmental clues (temperature and photo- 
period) that occur during an annual cycle. This process 
entails a fine modulation of a vast array of genes in order 
to provide allostasis. At a cellular level, one of the most 
dramatic differences corresponds to the reorganization 
of the nucleolar components observed during seasonal 
adaptation, when rRNA transcription and processing are 
clearly affected [45,46]. In this context, we previously 
reported that epigenetic mechanisms are clearly involved 
in the phenotypical plasticity that is achieved by the 
carp. In particular, the hypermethylated state of riboso- 
mal gene promoter in winter-acclimatized carp, as well 
as histone H2A replacement by variant subtypes of 
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Figure 5 Seasonal enrichment analysis of H2A.Z on carp ribosomal cistron. (A) Schematic representation of genomic organization of the 
carp ribosomal cistron. Green stars indicate the regions of the qPCR amplicons associated in ChIP experiments. (B) The seasonal enrichment of 
H2A.Z in rDNA cistron is shown. (C) The seasonal enrichment of H2A.Z on RNA polymerase II promoters {L41, A 9 -desoturose, and JAK1) is shown 
(n = 3 in each season). The data presents real-time PCR measurements of the immunoprecipitated DNA at the corresponding rDNA regions. 
Standard deviations (± SD) are shown (*P <0.001, **P <0.005, Student's f-test). (D) Seasonal enrichment of canonical histone H3 in rDNA cistron. 
ChIP, chromatin immunoprecipitation; qPCR, quantitative PCR; SD, standard deviation. 



macroH2A, is dependent on the seasonal transcriptional 
activity of ribosomal genes, corroborating that changes 
in the architecture of the chromatin are reversible and 
dynamically respond to environmental changes [32,33]. 

In the present study, the expression of H2A.Z was 
detected in all tissues studied. However, we observed 
that the total content of H2A.Z transcripts can change 
in the liver and testis in a seasonally- dependent manner. 
The same difference was found in western blot experi- 
ments, when liver H2A.Z content tended to be higher 
during the winter than summer. On the other hand, in 
immunodetection assays from pituitary cells, protein 
content of H2A.Z was significantly higher in winter- 
acclimatized fish. Nevertheless, this result is not corre- 
lated with mRNA levels. It is important to notice that in 
the present study, polyclonal antibodies were used for 
immunodetection and these antibodies cannot differenti- 
ate among different H2A.Z subtypes because they are di- 
rected to the extremely conserved C-terminal region of 
the protein; thus, limiting us to relate mRNA levels with 
protein levels of different H2A.Z subtypes. 

We report for the first time four distinct subtypes of 
histone H2A.Z in the same organism. One of the de- 
duced amino acid sequences from carp shares a 100% 
identity with the H2A.Z.2.1 of humans [24], whereas the 



other two (H2A.Z.3.1 and H2A.Z.3.2) do not correspond 
to any H2A.Z described until now. Between them, they 
share a 99.2% identity, differing only in Glyll to Ser. In 
carp, we could not identify the existence of histone 
H2A.Z.1, but we do not discard its existence in other tis- 
sues different to those studied in this work. 

Unexpectedly, we identified a new subtype of the his- 
tone variant H2A.Z (which was designated as H2A.Z.7) 
that is similar to H2A.Z.3.2, but possesses seven additional 
residues, more than other H2A.Z variants. The extra co- 
dons in H2A.Z.7 mRNA are in the same location as an in- 
tron in human H2A.Z, most likely indicating that H2A.Z7 
is a splicing (or mis-splicing) variant of H2A.Z.3.2. These 
seven additional residues are located in the region that 
corresponds to a-helix 2 [47], which could deviate H2A.Z 
C-terminal docking domain, most likely affecting the 
interaction surface with (H3-H4) 2 tetramer. Clearly, struc- 
tural studies should be carried out to determine the func- 
tional implications of these additional residues. However, 
considering we have no indication that the H2A.Z.7 pro- 
tein is translated, it could also be possible that this corre- 
sponds to a pseudogene. 

H2A.Z.7 presents a tissue-specific expression, which 
was detected only in the brain (summer and winter) and 
pituitary tissues (winter). This is not the first time that 
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Figure 6 Deposition of H2A.Zub on carp ribosomal cistron and RNA polymerase II promoters. Sequential chromatin immunoprecipitation 
(ChlP-reChIP) assays showing (A) seasonal enrichment of H2A.Zub along carp ribosomal cistron, and (B) seasonal enrichment of H2A.Zub on the 
L41 and A 9 -desaturase promoters. (C) Seasonal average ratio (summer/winter) of mRNA levels of ubiquitylation/deubiquitylation factors. The data 
were normalized with respects to carp 16S mitochondrial rRNA (n = 3 in each season). Standard deviations (± SD) are shown (*P = 0.047, 
**P = 0.002, ***P = 0.523, Student's f-test). ChlP-reChIP, sequential chromatin immunoprecipitation; SD, standard deviation. 



an unusual length of H2A.Z has been described In human 
tissues and cell lines, a shorter version named H2A.Z.2.2 
was found, and whose incorporation leads to severely un- 
stable nucleosomes due to it containing a characteristic- 
docking domain [24]. 

Compared to mammals, the presence of this particular 
and species-specific histone H2A.Z.7 could be due to 
the three whole genome duplications (WGDs) that tele- 
ost have experienced. In particular, carp are believed to 
have had another round of genome duplications (4R) 
and became an evolutionarily recent tetraploid fish [48]. 
This suggests that the presence of H2A.Z.7, and the 
other forms of carp H2A.Z, could be a consequence of 
WGD as a selective and important mechanism for 
eukaryote genome evolution in response to biotic and 
abiotic stress in carp. 

The clustering pattern of carp H2A.Z subtypes in the 
phylogenetic tree is consistent with the differentiations 
and diversifications of teleost fish. So far, H2A.Z.7 from 
carp has an intense evolutionary diversification between 
other H2A.Z from mammals and, even more, within 
teleost fish. The presence of more than one subtype of 
H2A.Z is not a whim of random genetic drift, but differ- 
ent forms may have acquired new or complementary 
functions [23]. Future experiments must show if the new 



subtypes described in this work perform similar or dif- 
ferent roles in chromatin. 

The incorporation of histone H2A.Z within chro- 
matin is crucial for proper gene regulation [17,22]. 
Genome-wide studies in several organisms show that 
H2A.Z is found throughout the genome; however, 
H2A.Z deposition is not random on chromatin, but ra- 
ther specific and highly regulated [35,49]. In yeast, 
H2A.Z is distributed in active and inactive RNA poly- 
merase II promoters [35,38,39], whereas in higher eu- 
karyotes it is mostly in active genes [50]. Interestingly, 
we report for the first time the seasonally-dependent 
enrichment of H2A.Z throughout the ribosomal cis- 
tron, concomitant with differential modulation of 
rDNA activity using a natural model of environmental 
plasticity, such as the carp fish. We observed that the 
levels of H2A.Z enrichment were higher in the riboso- 
mal cistron during summer, a season in which rRNA 
transcription is highly active [45]. This particular evi- 
dence confirms that the correlation between H2A.Z 
enrichment with the ribosomal transcriptional activity 
is not only associated to RNA polymerase I-related 
genes, since we detected that the enrichment of H2A.Z 
agrees with seasonally active gene promoters that are 
transcribed by RNA polymerase II. 
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Figure 7 Colocalization of H2A.Z with epigenetic markers during carp seasonal acclimatization in ribosomal cistron, and L41 and 
A9-desaturase gene promoters. Sequential chromatin immunoprecipitation (ChlP-reChIP) assays showing the seasonal colocalization of H2A.Z 
with (A) the euchromatin marker (H4K12ac) and (B) heterochromatin marker (H3K9me3) along the ribosomal cistron. (C,D) Represent the same 
seasonal colocalization, but in L41 and A 9 -desaturose promoters (n = 3 for each season). Standard deviations (± SD) are shown (*P <0.005, 
**P <0.001, ***P = 0.019, ****p = 0.015, Student's t-test). ChlP-reChIP, sequential chromatin immunoprecipitation; SD, standard deviation. 



In higher eukaryotes, p400 and SRCAP (SWI2/SNF2- 
related CBP activator protein) complexes are ATP- 
dependent chromatin remodelers specifically exchanging 
canonical H2A-H2B with H2A.Z-H2B dimers within 
the nucleosomes [51,52]. In carp, both glycogen and 
lipid content in hepatocytes from acclimatized carp dif- 
fer significantly between seasons [53]. Taking into ac- 
count that ribosome production is a major biosynthetic 
and energy-consuming activity of eukaryotic cells, which 
adapts rapidly to changes in intracellular energy status 
[54], it is plausible to speculate that variations of cellular 
AMP/ATP ratios could modulate H2A.Z-H2B incorpor- 
ation by p400 and SRCAP complexes in an ATP- 
availability manner. Certainly, this speculation must be 
demonstrated experimentally. 

Additionally, in order to corroborate if H2A.Z enrich- 
ment is consistent with a chromatin state, we evaluated 
H2A.Z colocalization with specific epigenetic markers of 
active (H4K12ac) [55] or inactive (H3K9me3) [56] gene 
expression. Our results showed that H2A.Z and the eu- 
chromatin marker (H4K12ac) co-localize along the ribo- 
somal cistron but not in a seasonally-dependent manner. 
In contrast, we showed a significant colocalization be- 
tween H3K9me3 and H2A.Z on rDNA cistron during 
winter, a season in which rRNA transcription and 



processing are highly inactive. Then, when we compared 
the colocalization of H2A.Z with these epigenetic markers 
on RNA polymerase II genes, they were consistent with 
the active or inactive gene states. These results agree with 
the controversial role of H2A.Z in the transcription of 
RNA polymerase II genes [57] . 

Histone tails are subjected to PTMs that alter the 
higher-order folding of chromatin. In particular, H2A.Z 
has been reported to be acetylated, sumoylated, and 
ubiquitylated. The H2A.Zub marker is involved in the 
maintenance or formation of heterochromatin because it 
was found predominantly on the inactive X chromosome 
[27]. In this study, we investigated if H2A.Z ubiquitylation 
in the ribosomal cistron can contribute to regulate carp 
seasonal rRNA transcription. First, we evaluated transcrip- 
tional expression of the most relevant ubiquitylation fac- 
tors. We found a greater gene expression of RING lb and 
Bmil in summer than in winter, whereas USP10 showed 
no change between seasons. These findings are consistent 
with the relevant enrichment of H2A.Zub on the riboso- 
mal cistron of summer- acclimatized carp, when the activ- 
ity of rRNA transcription is significantly higher. Unlike 
the uniform pattern of H2A.Z on the ribosomal cistron 
during summer, H2A.Zub enrichment in the same season 
was found mainly on the IGS and CP regions. These two 
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locations of the rDNA gene might play a crucial role in 
rRNA transcription since H2A ubiquitylation (and conse- 
quently H2A.Zub) can serve as a docking site for recruiting 
other factors, or preventing chromatin access to transcrip- 
tional regulators [58]. On the other hand, at least in hepa- 
tocytes isolated from carp, the two promoters of RNA 
polymerase II studied did not show significant enrichment 
of H2A.Zub in any season. Nevertheless, similar findings 
have been reported in Saccharomyces cerevisiae, where 
H2BK123ub has negative consequences on the assembly of 
RNA polymerase II at promoters and stimulates transcrip- 
tion elongation [59]. Finally, it would certainly be desirable 
to obtain antibodies specific for carp H2A.Zub to unam- 
biguously corroborate the ubiquitylation state of carp 
H2A.Z. Additionally, wide-genome screening of H2A.Zub 
should be performed in order to better understand the 
contribution of this PTM in transcriptional regulation dur- 
ing carp seasonal acclimatization. 

Conclusions 

Our findings reveal seasonal H2A.Z enrichment in active 
genes of the ribosomal cistron and RNA polymerase II 
promoters. In the same context, it seems that H2A.Zub 
enrichment is a specific marker of post-translational his- 
tone modification on active ribosomal cistron during the 
carp acclimatization process. Furthermore, we identified 
for the first time four distinct subtypes for the histone 
H2A.Z in the same organism. Further research is neces- 
sary to distinguish the particular contribution of each of 
the H2A.Z subtype within chromatin, and how the pre- 
dicted monoubiquitylation of H2A.Z modulates gene 
transcription. 

Methods 

Animal and tissue preparation 

All experimental work on fish described in this paper was 
approved by the Bioethics Committee of the Universidad 
Andres Bello, Santiago, Chile. Male carp were maintained 
in seasonal and environmental conditions as previously 
described [32]. The fish were sacrificed, and tissue was 
collected from the liver, brain, pituitary, testis, and gills, 
washed with PBS, pH 7.4, and stored at -80°C until 
required. 

Isolation of carp H2A.Z sequences 

From Danio rerio, H2A.Z mRNA [NCBI Reference Se- 
quence:NM_001201563.1] were deduced as the following 
oligonucleotides: Z-380pb, F 5 -GCAGGTGGAAAAGCA 
GGT AAAG-3 ' and Z-380pb, R 5 -GCGGTTTTCTGCT 
GGCCCTTCTT-3'. A 380 bp fragment corresponding to 
a partial mRNA region of H2A.Z was amplified by RT- 
PCR using carp liver tissue as a template. The product 
was cloned in the pGEM-T vector (Promega, Madison, 
WI, USA) and sequenced. Full mRNA of different carp 



H2A.Z subtypes were isolated using FirstChoice RLM- 
RACE (Ambion, Foster City, CA, USA) with specific 
oligonucleotides deduced from the fragments described 
earlier. The 5 'UTR region was completed with Z-380pb, R 
and 5'Z inner 5 -CCAGCTGTAGATGA CGAG-3', 
whereas the 3 'UTR region was completed with Z-380 pb, 
F and 3'Z inner 5 -GTCTACAGTGCAGCCATC-3'. All 
PCR reactions were performed with Phusion High-Fidelity 
DNA Polymerase (New England Biolabs, Ipswich, MA, 
USA) and DNA sequencing services were performed by 
Macrogen (Tokyo, Japan). 

Immunofluorescence 

A total of 100 mg of pituitary tissue from summer- and 
winter-acclimatized carp were fixed using 4% parafor- 
maldehyde (w/v). Sections were embedded in Paraplast 
and cut at 5 (im in a cryostat Leica CM1850 (Leica 
Microsystems, Wetzlar, Germany). The primary antibody 
used was anti-H2A.Z (1:500, ab4174; Abeam, Cambridge, 
MA, USA) incubated at 4°C for 16 hours and then washed 
three times with PBS. The secondary antibody was anti- 
rabbit IgG (1:400) conjugated to Alexa Fluor 647 (Cell 
Signaling Technology, Danvers, MA, USA) and 4',6- 
diamidino-2-phenylindole (DAPI) (1:1000), which was 
incubated for 1 hour at 25°C. Confocal microscopy was 
carried out in Axiovert 100 M (Zeiss, Oberkochen, 
Germany) and analyzed with ImageJ software (National 
Institutes of Health (NIH), Bethesda, MD, USA). 

RT-PCR and RT-qPCR 

Total RNA was purified from tissues, such as the liver, 
brain, pituitary, testis, and gills using TRIzol reagent 
(Invitrogen, Carlsbad, CA, USA). The quality of RNA was 
evaluated on 1.2% agarose gel electrophoresis and quanti- 
fied spectrophotometrically at A 2 60nm- The cDNA was 
synthesized from 2 |ig of total RNA with M-MLV Reverse 
Transcriptase (Promega) and random primers according 
to the manufacturers instructions. The resulting cDNAs 
were PCR-amplified in GoTaq Green Master Mix 
(Promega) using oligonucleotides, showed in Additional 
file 4: Table SI. 

RT-qPCR was performed with the appropriate oligo- 
nucleotides (Additional file 4: Table SI) and Brilliant 
SYBR Master Mix (Stratagene, La Jolla, CA, USA) in a 
final volume of 20 \i\ on a Mx3000P Real-Time PCR Sys- 
tem (Stratagene). Seasonal average ratio of mRNA levels 
was calculated using the Pfaffl method [60] and carp 16S 
mitochondrial rRNA as a seasonally constitutive gene 
expression. 

Western blot 

Nuclear extracts from HeLa cells, and rat and carp hepa- 
tocytes were prepared according to Araya et al. [33]. For 
the immune detections, we used polyclonal anti-H2A.Z 
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(1:1500; Abeam) and anti-H2B (1:5000; Abeam) antibodies. 
The peroxidase-conjugated anti-rabbit IgG (H+L) second- 
ary antibodies were obtained from KPL (Gaithersburg, 
MD, USA). Protein-antibody complexes were visualized by 
an enhanced chemiluminescence (ECL) detection system 
(GE Healthcare, Piscataway, NJ, USA). 

Chromatin immunoprecipitation combined with qPCR 

ChIP and sequential chromatin immunoprecipitation 
(ChlP-reChIP) assays were undertaken as previously de- 
scribed [33]. For each assay, we used 3 x 10 7 nuclei iso- 
lated from the liver of summer- and winter-acclimatized 
carp, respectively. ChlPs were performed with polyclonal 
anti-H2A.Z (ab4174; Abeam) and anti-H3 (abl791; 
Abeam). The ChlP-reChIP assays were performed immu- 
noprecipitating first with anti-H2A.Z antibodies, and then 
an additional immunoprecipitation with polyclonal his- 
tone H4 (acetyl K12), anti-histone H3 (Tri-methyl K9; 
Abeam), and anti-ubiquitin (Cell Signaling Technology). 
In the last case, although we did not have antibodies for 
carp H2A.Zub, our results were in agreement with H2A. 
Zub form because all data values obtained were higher 
than the unspecific IgG IP (normalizer). The immu- 
noprecipitated DNA was amplified by qPCR on an 
Mx3000P Real-Time PCR System (Stratagene) using oligo- 
nucleotides corresponding to IGS, TV , T 0 , CP and 5.8S re- 
gions of carp rDNA (Figure 5A), and the promoter 
regions of L41 [GenBank:AYl 15478], A 9 -desaturase 
[GenBank:BX571725.1] and JAK1 [GenBank:AH004872.1] 
genes (Additional file 4: Table SI). 

ChlP-qPCR was performed using SYBR Master Mix 
(Stratagene). The specificity of the PCR reaction was 
assessed by melting-temperature profiles. The fold en- 
richment was calculated using the AAQ method [61]: 

Fold Enrichment = 2^ AACt PUP/^) 
^[normalized Chip] = (Q X -(Q input] -Log 2 (input dilution factor))) 

[normalized IgG] = Q IgG "Log 2 (input dilution factor) 

^Q[ChIP/IgG] = [normalized ChIP] ~A Q [normalized IgG] 

Where C tx corresponds to the Ct value for each spe- 
cific antibody used. All assays were performed in tripli- 
cate and repeated with three different animal samples 
for each season. Results represent the mean (± SD) of 
three experiments. 

Phylogenetic analysis of H2A.Z variants 

The CDS from different eukaryotes corresponding to H2A. 
Z were retrieved from the GenBank database (Additional 
file 5: Table S2). A total of 29 sequences were aligned using 
ClustalW alignment [62] and analyzed by phylogeny recon- 
struction using the minimum evolution (ME) method and 
the Close-Neighbor-Interchange (CNI) algorithm. The tree 
is drawn to scale, with branch lengths in the same units as 
those of the evolutionary distances used to infer the 



phylogenetic tree. The evolutionary distances were com- 
puted using the p-distance method. The bootstrap consen- 
sus tree inferred from 5,000 replicates was taken to 
represent the evolutionary history of the taxa analyzed. 
Evolutionary analyses were conducted in MEGA 5.05 soft- 
ware (Tokyo, Japan) [63]. 

Additional files 



Additional file 1: Figure SI. Carp fish exposed to a wide range of 
temperature and photoperiod conditions. Graph of temperature and 
photoperiod changes during a seasonal annual cycle. 

Additional file 2: Figure S2. Identification and analysis of H2A.Z 
subtypes isolated from carp. Multiple nucleotide and amino acid 
sequence alignment of H2A.Z subtypes isolated from carp are shown. 
Black boxes represent the precise sites of nucleotide differences between 
CDS from H2A.Z.2.1, H2AZ3.1, H2A.Z.3.2, and H2A.Z.7. The deduced 
amino acid sequences of each H2A.Z subtype are shown in colors. The 
additional seven amino acids found in carp H2A.Z.7 are denoted with a 
red rectangle. CDS, coding DNA sequences. 

Additional file 3: Figure S3. Constitutive gene expression of JAK1 
during seasonal acclimatization of carp. Transcriptional expression analysis 
of JAK1 in liver tissues from summer- and winter-acclimatized carp. The 
mRNA levels of JAK1 were calculated using the Pfaffl method [60]. Three 
different fish for each season were used. 

Additional file 4: Table SI. List of oligonucleotides used for RT-PCR, 
RT-qPCR, and ChIP assays. ChIP, chromatin immunoprecipitation; RT-PCR, 
reverse transcription polymerase chain reaction; RT-qPCR, reverse 
transcription quantitative PCR. 

Additional file 5: Table S2. GenBank Accession Numbers of CDS used 
for the phylogenetic analysis. The list depicts histone H2A.Z in eukaryotic 
organisms used in the phylogenetic analysis. CDS, coding DNA 
sequences. 



Abbreviations 

bp: base pairs; CDS: Coding DNA sequences; ChIP: Chromatin 
immunoprecipitation, ChlP-reChIP, sequential chromatin 
immunoprecipitation; CNI: Close-Neighbor-Interchange; CP: Core promoter; 
DAPI: 4',6-diamidino-2-phenylindole; ECL: Enhanced chemiluminescence; 
H2A.Zub: ubiquitylated-H2A.Z; H3K9me3: Trimethylated histone H3 at lysine 
9; H4K12ac: acetylated histone H4 at lysine 12; IgG: Immunoglobulin G; 
IGS: Intergenic spacer; ME: Minimum evolution; PBS: Phosphate buffered 
saline; PCR: Polymerase chain reaction; PTMs: Post-translational modifications; 
qPCR: quantitative PCR; RT: Reverse transcription; RT-PCR: Reverse 
transcription polymerase chain reaction; SD: Standard deviation; 
UCD: University College Dublin; WGD: Whole genome duplication. 

Competing interests 

The authors declare that they have no competing interests. 
Authors' contributions 

NS carried out the experiments and wrote the manuscript. GN carried out 
some ChIP assays. MR carried out the immunofluorescence experiment. 
* MR is currently working at Departamento Ciencias Fundantes, Universidad 
Central de Chile. MA directed the research and co-wrote the manuscript. 
AM assisted with project design and development, and drafted the 
manuscript. All authors read and approved the final manuscript. 

Acknowledgements 

Financial support was received from Fondo Nacional de Desarrollo Cientifico 
yTecnologico (FONDECY~0: 1120873; Universidad Andres Bello (UNAB): 
DI-451 1/R. MA also acknowledges to the Comision Nacional de Investigacion 
Cientifica y Tecnologica (CONICYT)/Fondo de Financiamiento de Centros de 
Investigacion en Areas Prioritarias (FONDAP): 15110027. We are grateful to 
Eduardo Fuentes for constructive comments on the manuscript, and to 
P Rudolffi and P Carvajal for technical support. 



Simonet et al. Epigenetics & Chromatin 2013, 6:22 
http://www.epigeneticsandchromatin.eom/content/6/1/22 



Page 12 of 13 



Author details 

1 Laboratories de Biologia Celular y Molecular, Facultad de Ciencias Biologicas, 
Universidad Andres Bello, Los Fresnos 52, Vina del Mar, Chile, 
interdisciplinary Center for Aquaculture Research (INCAR), Victor Lamas 
1290, P.O. Box 160-C, Concepcion, Chile. 

Received: 28 November 2012 Accepted: 7 June 2013 
Published: 17 July 2013 

References 

1. Lopez-Maury L, Marguerat S, Bahler J: Tuning gene expression to changing 
environments: from rapid responses to evolutionary adaptation. 

Nat Rev Genet 2008, 9:583-593. 

2. Cossins AR, Crawford DL: Fish as models for environmental genomics. Nat 
Rev Genet 2005, 6:324-333. 

3. Alvarez M, Molina A, Quezada C, Pinto R, Krauskopf M, Vera M: Eurythermal 
fish acclimatization and nucleolar function: a review. J Therm Biol 2004, 
29:663-667. 

4. Vaquero A, Loyola A, Reinberg D: The constantly changing face of 
chromatin. Sci Aging Knowledge Environ 2003, 2003:RE4. 

5. Luger K, Mader AW, Richmond RK, Sargent DF, Richmond TJ: Crystal 
structure of the nucleosome core particle at 2.8 A resolution. Nature 
1997, 389:251-260. 

6. Bollati V, Baccarelli A: Environmental epigenetics. Heredity 201 0, 1 05:1 05-1 1 2. 

7. Rottach A, Leonhardt H, Spada F: DNA methylation-mediated epigenetic 
control. J Cell Biochem 2009, 108:43-51. 

8. Peterson CL: ATP-dependent chromatin remodeling: going mobile. 
FEB5 Lett 2000, 476:68-72. 

9. Esteller M: Non-coding RNAs in human disease. Nat Rev Genet 201 1, 
12:861-874. 

10. Henikoff S, Ahmad K: Assembly of variant histones into chromatin. 

Annu Rev Cell Dev Biol 2005, 21 :1 33-1 53. 

1 1 . Kouzarides T: Chromatin modifications and their function. Cell 2007, 
128:693-705. 

12. Jin J, Cai Y, Li B, Conaway RC, Workman JL, Conaway JW, Kusch T: In and 
out: histone variant exchange in chromatin. Trends Biochem Sci 2005, 
30:680-687. 

13. Kamakaka RT: Histone variants: deviants? Genes Dev 2005, 19:295-316. 

14. louzalen N, Moreau J, Mechali M: H2A.ZI, a new variant histone expressed 
during Xenopus early development exhibits several distinct features 
from the core histone H2A. Nucleic Acids Res 1 996, 24:3947-3952. 

15. Thatcher TH, Gorovsky MA: Phylogenetic analysis of the core histones 
H2A, H2B, H3, and H4. Nucleic Acids Res 1994, 22:174-179. 

16. Jackson J, Gorovsky MA: Histone H2A.Z has a conserved function that is 
distinct from that of the major H2A sequence variants. Nucleic Acids Res 
2000, 28:3811-3816. 

17. Clarkson MJ, Wells JR, Gibson F, Saint R, Tremethick DJ: Regions of variant 
histone His2AvD required for Drosophila development. Nature 1999, 
399:694-697. 

18. Liu X, Bowen J, Gorovsky MA: Either of the major H2A genes but not an 
evolutionarily conserved H2A.F/Z variant of Tetrahymenathermophila 
can function as the sole H2A gene in the yeast Saccharomyces 
cerevisiae. Mol Cell Biol 1996, 16:2878-2887. 

19. Ridgway P, Brown KD, Rangasamy D, Svensson U, Tremethick DJ: Unique 
residues on the H2A.Z containing nucleosome surface are important for 
Xenopuslaevis development. J Biol Chem 2004, 279:43815-43820. 

20. Faast R, Thonglairoam V, Schulz TC, Beall J, Wells JR, Taylor H, Matthaei K, 
Rathjen PD, Tremethick DJ, Lyons I: Histone variant H2A.Z is required for 
early mammalian development. Curr Biol 2001, 1 1:1 183-1 187. 

21. Nashun B, Yukawa M, Liu H, Akiyama T, Aoki F: Changes in the nuclear 
deposition of histone H2A variants during pre-implantation 
development in mice. Development 2010, 137:3785-3794. 

22. Adam M, Robert F, Larochelle M, Gaudreau L: H2A.Z is required for global 
chromatin integrity and for recruitment of RNA polymerase II under 
specific conditions. Mol Cell Biol 2001, 21:6270-6279. 

23. Einn-Lopez JM, Gonzalez-Romero R, Dryhurst D, Ishibashi T, Ausio J: The 
evolutionary differentiation of two histone H2A.Z variants in chordates 
(H2A.Z-1 and H2A.Z-2) is mediated by a stepwise mutation process that 
affects three amino acid residues. BMC Evol Biol 2009, 9:31. 

24. Bonisch C, Schneider K, Punzeler S, Wiedemann SM, Bielmeier C, Bocola M, 
Eberl HC, Kuegel W, Neumann J, Kremmer E, Leonhardt H, Mann M, 



Michaelis J, Schermelleh L, Hake SB: H2A.Z.2.2 is an alternatively spliced 
histone H2A.Z variant that causes severe nucleosome destabilization. 

Nucleic Acids Res 201 2, 40:5951-5964. 

25. Millar CB, Xu F, Zhang K, Grunstein M: Acetylation of H2AZ Lys 14 is 
associated with genome-wide gene activity in yeast. Genes Dev 2006, 

20:711-722. 

26. Kalocsay M, Hiller NJ, Jentsch S: Chromosome-wide Rad51 spreading and 
SUMO-H2A.Z-dependent chromosome fixation in response to a 
persistent DNA double-strand break. Mol Cell 2009, 33:335-343. 

27. Sarcinella E, Zuzarte PC, Lau PN, Draker R, Cheung P: Monoubiquitylation of 
H2A.Z distinguishes its association with euchromatin or facultative 
heterochromatin. Mol Cell Biol 2007, 27:6457-6468. 

28. Ben-Saadon R, Zaaroor D, Ziv T, Ciechanover A: The polycomb protein 
Ringl B generates self atypical mixed ubiquitin chains required for its 
in vitro histone H2A ligase activity. Mol Cell 2006, 24:701-71 1. 

29. Cao R, Tsukada Y, Zhang Y: Role of Bmi-1 and Ringl A in H2A 
ubiquitylation and Hox gene silencing. Mol Cell 2005, 20:845-854. 

30. Draker R, Sarcinella E, Cheung P: USP10 deubiquitylates the histone 
variant H2A.Z and both are required for androgen receptor-mediated 
gene activation. Nucleic Acids Res 201 1, 39:3529-3542. 

31. Grummt I, Pikaard CS: Epigenetic silencing of RNA polymerase I 
transcription. Nat Rev Mol Cell Biol 2003, 4:641-649. 

32. Pinto R, Ivaldi C, Reyes M, Doyen C, Mietton F, Mongelard F, Alvarez M, 
Molina A, Dimitrov S, Krauskopf M, Vera Ml, Bouvet P: Seasonal 
environmental changes regulate the expression of the histone 
variant macroH2A in an eurythermal fish. FEBS Lett 2005, 
579:5553-5558. 

33. Araya I, Nardocci G, Morales J, Vera M, Molina A, Alvarez M: MacroH2A 
subtypes contribute antagonistically to the transcriptional regulation of 
the ribosomal cistron during seasonal acclimatization of the carp fish. 

Epigenetics Chromatin 2010, 3:14. 

34. Vera Ml, Molina A, Pinto R, Reyes M, Alvarez M, Krauskopf E, Quezada C, 
Torres J, Krauskopf M: Genomic organization of the rDNAcistron of the 
teleost fish Cyprinuscarpio. Biol Res 2003, 36:241-251. 

35. Zhang H, Roberts DN, Cairns BR: Genome-wide dynamics of Htz1, a 
histone H2A variant that poises repressed/basal promoters for activation 
through histone loss. Cell 2005, 123:219-231. 

36. Barski A, Cuddapah S, Cui K, Roh T, Schones DE, Wang Z, Wei G, Chepelev I, 
Zhao K: High-resolution profiling of histone methylations in the human 
genome. Cell 2007, 129:823-837. 

37. Mavrich TN, Jiang C, loshikhes IP, Li X, Venters BJ, Zanton SJ, Tomsho LP, Qi 
J, Glaser RL, Schuster SC, Gilmour DS, Albert I, Pugh BF: Nucleosome 
organization in the Drosophila genome. Nature 2008, 453:358-362. 

38. Guillemette B, Bataille AR, Gevry N, Adam M, Blanchette M, Robert FO, Gaudreau 
L: Variant histone H2A.Z is globally localized to the promoters of inactive 
yeast genes and regulates nucleosome positioning. PLoS Biol 2005, 3:e384. 

39. Li B, Pattenden SG, Lee D, Gutierrez J, Chen J, Seidel C, Gerton J, Workman 
JL: Preferential occupancy of histone variant H2AZ at inactive promoters 
influences local histone modifications and chromatin remodeling. 

Proc Natl Acad Sci USA 2005, 1 02:1 8385-1 8390. 

40. Hardy S, Jacques P, Gevry N, Forest A, Fortin M, Laflamme L, Gaudreau L, 
Robert F: The Euchromatic and heterochromatic landscapes are shaped 
by antagonizing effects of transcription on H2A.Z Deposition. PLoS Genet 
2009, 5:e1 000687. 

41. Valdes-Mora F, Song JZ, Statham AL, Strbenac D, Robinson MD, Nair SS, 
Patterson Kl, Tremethick DJ, Stirzaker C, Clark SJ: Acetylation of H2A.Z is a 
key epigenetic modification associated with gene deregulation and 
epigenetic remodeling in cancer. Genome Res 2012, 22:307-321. 

42. Talbert PB, Ahmad K, Almouzni G, Ausio J, Berger F, Bhalla PL, Bonner WM, 
Cande WZ, Chadwick BP, Chan SWL, Cross GAM, Cui L, Dimitrov SI, 
Doenecke D, Einn-Lopez JM, Gorovsky MA, Hake SB, Hamkalo BA, Holec S, 
Jacobsen SE, Kamieniarz K, Khochbin S, Ladurner AG, Landsman D, Latham 
JA, Loppin B, Malik HS, Marzluff WF, Pehrson JR, Postberg J, et al: A unified 
phylogeny-based nomenclature for histone variants. Epigenetics Chromatin 
2012,5:7. 

43. Gracey A, Macartney A, Beynon R, Cossins A: Cold-induced expression of 
A 9 -desaturase in carp by transcriptional and posttranslational 
mechanisms. Science 1996, 271:815-818. 

44. Molina A, Corta A, Martin R, Alvarez M, Burzio LO, Krauskopf M, Vera Ml: 
Gene structure of the carp fish ribosomal protein L41: seasonally 
regulated expression. Biochem Biophys Res Commun 2002, 295:582-586. 



Simonet et al. Epigenetics & Chromatin 2013, 6:22 
http://www.epigeneticsandchromatin.eom/content/6/1/22 



Page 13 of 13 



45. 



46. 



47. 



50. 



52. 



53. 



54. 



55. 



56. 



57. 



58. 



59. 



60. 



62. 



63. 



Vera Ml, Norambuena L, Alvarez M, Figueroa J, Molina A, Leon G, Krauskopf 
M: Reprog ramming of nucleolar gene expression during the 
acclimatization of the carp. Cell Mol Biol Res 1993, 39:665-674. 
Alvarez M, Quezada C, Molina A, Krauskopf M, Vera Ml, Thiry M: 
Ultrastructural changes of the carp (Cyprinus carpio) hepatocyte 
nucleolus during seasonal acclimatization. Biol Cell 2006, 98:457-463. 
Suto RK, Clarkson MJ, Tremethick DJ, Luger K: Crystal structure of a 
nucleosome core particle containing the variant histone H2A.Z. 
Nat Struct Biol 2000, 7:1 1 21 -1 1 24. 

Danzmann RG, Davidson EA, Ferguson MM, Gharbi K, Koop BF, Hoyheim B, 
Lien S, Lubieniecki KP, Moghadam HK, Park J, Phillips RB, Davidson WS: 
Distribution of ancestral proto-Actinopterygian chromosome arms within 
the genomes of 4R-derivative salmonid fishes (Rainbow trout and 
Atlantic salmon). BMC Genomics 2008, 9:557. 
Billon P, Cote J: Precise deposition of histone H2A.Z in chromatin for 
genome expression and maintenance. Biochim Biophys Acta 2012, 
1819:290-302. 

Weber CM, Henikoff JG, Henikoff S: H2A.Z nucleosomes enriched over 
active genes are homotypic. Nat Struct Mol Biol 2010, 17:1500-1507. 
Ruhl DD, Jin J, Cai Y, Swanson S, Florens L, Washburn MP, Conaway RC, 
Conaway JW, Chrivia JC: Purification of a human SRCAP complex that 
remodels chromatin by incorporating the Histone variant H2A.Z into 
nucleosomes. Biochemistry 2006, 45:5671-5677. 
Cai Y: The mammalian YL1 protein is a shared subunit of the TRRAP/ 
TIP60 histone acetyltransferase and SRCAP complexes. J Biol Chem 2005, 
280:13665-13670. 

Saez L, Zuvic T, Amthauer R, Rodriguez E, Krauskopf M: Fish liver protein 
synthesis during cold acclimatization: seasonal changes of the 
ultrastructure of the carp hepatocyte. J Exp Zool 1 984, 230:1 75-1 86. 
Murayama A, Ohmori K, Fujimura A, Minami H, Yasuzawa-Tanaka K, Kuroda 
T, Oie S, Daitoku H, Okuwaki M, Nagata K, Fukamizu A, Kimura K, Shimizu T, 
Yanagisawa J: Epigenetic control of rDNA loci in response to intracellular 
energy status. Cell 2008, 133:627-639. 

Bernstein BE, Kamal M, Lindblad-Toh K, Bekiranov S, Bailey DK, Huebert DJ, 
McMahon S, Karlsson EK, Kulbokas EJ, Gingeras TR, Schreiber SL, Lander ES: 
Genomic maps and comparative analysis of histone modifications in 
human and mouse. Cell 2005, 120:169-181. 

Berletch JB, Yang F, Disteche CM: Escape from X inactivation in mice and 

humans. Genome Biol 2010, 11:213. 

Marques M, Laflamme L, Gervais AL, Gaudreau L: Reconciling the positive 
and negative roles of histone H2A.Z in gene transcription. Epigenetics 
2010,5:267-272. 

Braun S, Madhani HD: Shaping the landscape: mechanistic consequences 
of ubiquitin modification of chromatin. EMBO Rep 2012, 13:619-630. 
Fleming AB, Kao CF, Hillyer C, Pikaart M, Osley MA: H2B ubiquitylation 
plays a role in nucleosome dynamics during transcription elongation. 

Mol Cell 2008, 31:57-66. 

Pfaffl MW: A new mathematical model for relative quantification in real- 
time RT-PCR. Nucleic Acids Res 2001, 29:45e-45. 

Livak KJ, Schmittgen TD: Analysis of relative gene expression data using 
real-time quantitative PCR and the 2(-Delta Delta CCD) Method. 

Methods 2001, 25:402-408. 

Thompson JD, Higgins DG, Gibson TJ: CLUSTAL W: improving the 
sensitivity of progressive multiple sequence alignment through 
sequence weighting, position-specific gap penalties and weight matrix 
choice. Nucleic Acids Res 1 994, 22:4673-4680. 
Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: 
molecular evolutionary genetics analysis using maximum likelihood, 
evolutionary distance, and maximum parsimony methods. 
Mol Biol Evol 201 1 , 28:2731-2739. 



doi:1 0.1 186/1 756-8935-6-22 

Cite this article as: Simonet et al.: Epigenetic regulation of the ribosomal 
cistron seasonally modulates enrichment of H2A.Z and H2A.Zub in 
response to different environmental inputs in carp {Cyprinus carpio). 

Epigenetics & Chromatin 2013 6:22. 



Submit your next manuscript to BioMed Central 
and take full advantage of: 

• Convenient online submission 

• Thorough peer review 

• No space constraints or color figure charges 

• Immediate publication on acceptance 

• Inclusion in PubMed, CAS, Scopus and Google Scholar 

• Research which is freely available for redistribution 



Submit your manuscript at 
www.biomedcentral.com/submit 



o 



BioMed Central 



